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In this talk we present recent next-to-leading order results relevant for LHC phenomenology obtained with the 
GOLEM method. After reviewing the status of this Feynman diagrammatic approach for multi-leg one-loop calcu- 
lations we discuss three applications: the loop-induced process gg — » Z* Z* and the virtual corrections to the five 
and six point processes qq —> ZZg and uu -^ sscc. We demonstrate that our method leads to representations 
of such amplitudes which allow for efficient phase space integration. In this context we propose a reweighting 
technique of the leading order unweighted events by local K-factors. 



1. Introduction 

The Large Hadron Collider will explore our un- 
derstanding of fundamental interactions in the 
multi-TeV range. Apart from testing the Higgs 
mechanism which is the final cornerstone of the 
Standard Model (SM), also a plethora of exten- 
sions of the SM will be put under scrutiny. What- 
ever final states will be detected, the initial state 
will always consist of QCD partons. The pertur- 
bative description of such processes is necessar- 
ily plagued by renormalisation and factorisation 
scale uncertainties. Only if next-to-leading order 
(NLO) corrections are included the logarithmic 
dependence on these scales is tamed and one ar- 
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rives at sufficiently reliable predictions for various 
signal and background processes. For a discussion 
of what remains to be done see [1] . 

For a full next-to- leading order evaluation of an 
N-point process one has to combine virtual cor- 
rections with real-emission corrections using some 
infrared subtraction method. The tree-like pro- 
cesses can be evaluated by using standard leading 
order tools. Meanwhile, automated ways to deal 
with the IR subtractions based on the Catani- 
Seymour dipole approach [5j are also on the mar- 
ket [314] . The evaluation of the one-loop contribu- 
tion for N-point processes is not yet at this level 
of automation, although new ideas emerged and a 
lot of progress has been made recently in various 
directions |5l6l7l8l9llOllirT^ . 
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2. The GOLEM method 

The aim of our collaboration is to provide a 
tool which allows for a numerically stable evalua- 
tion of multi-leg one-loop amplitudes: GOLEPlj. It 
is based on the method described in ^ 13I14| . The 
approach relies on the evaluation of Feynman dia- 
grams and the reduction of tensor integrals using 
a form factor approach. The form factors can be 
evaluated in various ways as outlined below. 

We organize the evaluation of a one-loop am- 
plitude as follows: 

• generate Feynman diagrams using QGRAF 
[I5] or FeynArts 3.2 [17^. 

• separate and perform colour algebra 

• project on helicity amplitudes 

At this point two independent set-ups exist. 
Firstly, a completely symbolic reduction to stan- 
dard scalar integrals with up to four external legs 
can be obtained using FORM [TB] and MAPLE. 

+ Cbubl2 + Ctadli + Ti- ■ 



The respective coefficients are rational polynomi- 
als in Mandelstam variables. The extraction of 
the rational part, TZ, of the amplitude can be done 
separately [18]. As long as no efficient tools for 
the manipulation of multivariate rational poly- 
nomials are available, interactive user input is 
needed to produce sufficiently compact amplitude 
expressions in the purely symbolic approach. Sec- 
ondly, apart from the symbolic approach, we pro- 
vide a numerical tensor reduction. Schematically 
the amplitude is expressed in terms of form fac- 
tors which resemble Feynman parameter integrals 
with Feynman parameters in the numerator. 
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These form factors are implemented in a 
F0RTAN90 code and can be evaluated by numerical 
reduction and also by using one-dimensional in- 
tegral representations. The form factors were de- 
signed to avoid the occurrence of so-called Gram 



determinants which usually hamper a numerically 
stable evaluation of large Feynman diagrammatic 
expressions. Our method has been successfully 
applied to various calculations with up to six 
point functions ^19|20 21 22 . 

3. Applications for LHC phenomenology 

We discuss now recent evaluations of three loop 
amplitudes which are relevant in the context of 
Higgs searches at the LHC. 

3.1. The process gg -^ Z*Z* 

In |23|24|25] it has been shown that the gluon 
induced production of charged vector boson pairs 
accounts for about 30 percent of the background 
to Higgs searches in that channel after cuts. Al- 
though a similar calculation for neutral vector 
bosons has been performed a long time ago [26] . 
no public code is available which motivated us to 
redo this calculation using our method. 

As only 4-point functions are present, a sym- 
bolic expression for the amplitude could be ob- 
tained, where numerically dangerous denomina- 
tors have been cancelled algebraically. The ex- 
pressions are implemented in a flexible computer 
program GG2ZZ ^7 which also contains the pho- 
tonic contributions. The size of the gluon contri- 
bution to the ZZ cross section in relation to the 
quark induced part is included in the following 
table: 





a{pp -^ Z*ij*)Z*{-f*) -> iU'i') [fb] 
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GDLEM=General One- Loop Evaluator for Matrix elements. 



For the numerical results we use the follow- 
ing set of input parameters: Mw = 80.419 GeV, 
Mz = 91.188 GeV, Gf = 1.16639 10"^ GeV-^, 
Tz = 2.44 GeV. The electromagnetic coupling 
is defined in the G^ scheme. The pp cross sec- 
tions are calculated at yi = 14 TeV employing 
the CTEQ6L1 and CTEQ6M parton distribution 
functions at tree- and loop-level, for more de- 
tails see [28]. Applying standard cuts: pre > 20 
GeV,75 GeV< Mi+i- < 105 GeV, \rie\ < 2.5, 75 
GeV, we find that the gluon contribution ac- 
counts for 14% to the total pp -^ ZZ process. 
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The qq contribution was evaluated using MCFM 

m- 

The effect of the photon contribution can be 
seen best in the invariant mass distribution of the 
4 leptons in Fig [TJ Between the one- and two-Z 
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Figure 1. The invariant mass distribution of the 
4 leptons including photonic contributions. Only 
a minimal cut M^+g- > 5 GeV is applied. 



threshold the interference effects are sizable. 

3.2. The process PP -^ ZZ+jet 

This process is of relevance in the context of 
Higgs searches in the Higgs plus jet channel. We 
have obtained analytic expressions for all 36 par- 
tonic one-loop helicity amplitudes qqZZg — > 
which contribute to this process. The colour 
structure is simple, one finds three different colour 
structures. By applying projection operators to 
each Feynman diagram, reducible scalar prod- 
ucts between the loop momentum and external 
momenta can be expressed by inverse propaga- 
tors and cancelled. In this way only rank one 
five-point functions remain, together with two-, 
three- and four-point tensor integrals. We use 
the 'tHooft-Veltman scheme which needs an ac- 
companying prescription for 75. By splitting the 
7-algebra and the loop momentum in a 4- and (n- 
4)-dimensional part we have to add a finite coun- 
terterm ^ (1 — UsCp /i^) to the axial coupling to 



guarantee that the axial and vector part of the 
vector boson rcnormalise in the same way |30j . 
After UV renormalization, only IR poles remain. 
A finite expression can be obtained by adding the 
Catani-Seymour insertion operator, /(e) [2], to 
the result. We have integrated the resulting ex- 
pression for the LHC energy over the phase space 
using the cut pTjet > 100 GeV and find: 



ctlo 

O'LO+virt. 



1003.1 ± 0.4 fb 
899.0 ± 4.7 fb 



Here Mw = 80.403 GeV, Mz = 91.1876 GeV, 
Gf = 1.16637- 10^5 GeV^2 ^re used, a is evalu- 
ated in the G^ scheme and we used fiR — fip = 
Mz- For the LO result we use the CTEQ6L1 and 
for NLO the CTEgSM parton distribution func- 
tions, as is evaluated using the LHAPDF routines. 
Note that the number does not include the closed 
fermion loop contribution to this process as it 
turned out to be numerically irrelevant. 

We compare the scale dependence of the LO 
term and the NLO virtual part in Fig [21 The fig- 
ures indicate that the inclusion of the virtual cor- 
rections indeed stabilise the prediction for scales 
around 2Mz. Note that the real emission part 
of the NLO correction will add another iip de- 
pendent term which stems from the initial state 
singularities. As new colour channels are present 
at NLO one expects actually a deterioration of 
the scale dependence, as has been observed in the 
PP ^ W+W-+iet case [S]. 

Results for this related process have been pre- 
sented by two other groups already [31|32j . We 
have also evaluated this process and compared 
our evaluation with both groups for single phase 
space points. Perfect agreement was found [T]. 

3.3. The process uu -^ ssbb 

Signatures of beyond SM processes contain typ- 
ically leptons, jets and missing energy and one 
easily reaches large numbers of final state par- 
tons at the LHC. Up to now no complete 6-point 
process, which is related to a 2 ^- 4 kinematics, 
is evaluated at next-to-leading order in as- Note 
that progress in that direction has been reported 
for the process pp -^ bbti at this workshop |33|34j . 
An example for another relevant 6-point process 
emerges in two Higgs doublet models, which lead 
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Figure 2. Scale dependence of the LO and virtual 
next-to-leading order corrections. The two curves 
in each plot show diagonal and anti-diagonal vari- 
ation of the scales. 



in certain parameter regions predominantly to 4 
b-jets in the final state. To understand the related 
background in detail, the process PP — > bbbb has 
to be known at NLO. Note that at the LHC one 
can safely neglect the bottom mass, if realistic pr 
and b separation cuts are applied. The partonic 
amplitudes uu -^ ssbb and gg -^ ssbb are sufh- 
cient to predict this cross section in the massless 
limit. We report here on the successful evalua- 
tion of the virtual part of the first of these two 
amplitudes. 

In detail, the uu — > ssbb amplitude can be 
written in terms of six colour structures. Two 
independent helicity amplitudes, ^++++++ and 
^+++^ , are needed. They were evaluated in 



two completely independent ways. In one evalua- 
tion a fully symbolic reduction to scalar integrals 
was performed, in the other one each Feynman 
diagram was mapped to a form factor representa- 
tion and translated into a FQRTRAN90 code. Both 
calculations are highly automated such that the 
evaluation of other processes only needs the re- 
spective Feynman diagrammatic input and the 
specification of colour and helicity projections. In 
the given case one has to evaluate 25 pentagon 
and 8 hexagon diagrams. After UV renormal- 
ization and adding the IR insertion operator all 
poles in l/(n — 4) cancel and one is left with 
a finite expression which can be evaluated nu- 
merically. The F0RTRAN90 code was organized 
such that the reevaluation of algebraic terms was 
avoided by recursive organization of the expres- 
sions and caching. The evaluation time for one 
phase space point of the full amplitude, summed 
over helicities and colour is about 0.8 seconds on 
a 3.2 GHz Intel Pentium 4 processor. As the in- 
tegration over phase space can be trivially par- 
allelised this is sufficiently fast in what concerns 
the evaluation of distributions. 

As the evaluation of the amplitude needs a 
large number of numerical operations, one typ- 
ically observes numerical problems in parts of the 
phase space where denominators become small 
and form factors, respectively scalar integrals, are 
not linearly independent anymore. If one inte- 
grates directly the LO plus finite virtual correc- 
tions over the phase space, adaptive numerical 
integrators tend to sample phase space points in 
these critical phase space regions. This happens if 
the induced variations influence the result at the 
order of the accuracy goal. To avoid this kind 
of destabilisation we have applied the following 
method for integrating the virtual NLO correc- 
tions. 

We first evaluate the LO contribution over the 
target phase space 



f^LO 



dxfo{x) 



and generate unweighted events -&,=!, ...,jv- The 
latter are related to a parameter transformation 
a; — *■ y on phase space such that the new variables 
have constant density dy ^ dxfo{x)/aLO- Any 
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observable O can be estimated by distributing the 
unweighted events, Ej, into the respective bins. 

N f 

An estimate of the LO plus virtual corrections 
can now be obtained using the same set of events. 
The relation 

Oho+virtuai = / dxfi{x) = a LO / dyK{y) , 



The histograms in Fig. [3] are filled with 200,000 



where K — /1//0 is a local K-factor, implies 

N 

<yLo 



(O) 



LO-\-virt. 



--wT.xiE,)^m, 



which is a simple reweighting of a LO event sam- 
ple. For this purpose the LO events should be 
evaluated with NLO pdfs. In this way no integra- 
tion over the finite virtual corrections is needed, 
one simply has to evaluate the virtual corrections 
for each unweighted event which belongs to a 
specified observable. Of course it still can hap- 
pen that this evaluation is plagued by numeri- 
cal problems but it does not negatively affect the 
sampling of test points in integration methods. 
This method leads to a good estimate if the vir- 
tual corrections are sufficiently close to the LO 
distribution, such that the unweighted events are 
also representative for the LO-f- virtual differential 
cross section. Note that this has to be fulfilled 
for any observable for which perturbation theory 
is meaningful in the first place. 

To illustrate this method we show in Fig. [3] the 
effect of the virtual contribution on the distri- 
bution of the leading jet, i.e. the jet with the 
highest energy. To define the LO 4-jet observ- 
able we use the cuts: ptj > 50 GeV, \r]j\ < 3 
and AR = ^ Ac^P- -)- Ar?2 > 0.3. The LO cross 
section and the corresponding unweighted events 
were evaluated using WHIZARD [35] with CTEQ6M 
pdfs and the scale choice X]7=iPtj/4, M = 100 
GeV and one-loop running for Ug. The LO + 
finite virtual contribution was evaluated as de- 
scribed above. For the LO and LO-|-virtual con- 
tribution we obtain 



olo 

^ LO+virtual 



88.5 ± 0.2 fb 
69.0 ± 0.2 fb 
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Figure 3. The pr and rapidity distribution of the 
leading jet. The full line is the LO, the dashed line 
is obtained by adding the UV/IR finite contribu- 
tion of the virtual part of the NLO prediction, as 
defined in the text. 



unweighted events. When evaluating the local 
K-factors, less than 1% of all points showed an 
indication of numerical instability. These criti- 
cal points where simply reevaluated by using the 
quadruple precision version of our code. A simi- 
lar evaluation of the gg — > ssbb amplitude and the 
combination with the real emission corrections is 
in progress. 
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4. Conclusion 

In this talk we have presented recent results 
of the GOLEM collaboration. The implementa- 
tion of our method to evaluate Feynman dia- 
grammatic representations of amplitudes in sym- 
bolic/numerical computer programs has been 
completed in the context of one-loop amplitude 
evaluations relevant for the LHC. Here we pre- 
sented results for the process gg -^ Z*Z*, and 
the virtual corrections to PP —f ZZj and uu —>■ 
ssbb. We have proposed a new indirect inte- 
gration method for virtual corrections which is 
based on the evaluation of local K-factors for un- 
weighted events defined by the LO cross section. 
We conclude that our method is numerically ef- 
ficient and can provide predictions for multi-leg 
one loop processes at TeV colHders. 
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